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Leakage  current  is  one  of  the  critical  aspects  to  consider  for  overhead 
transmission  line  insulator’s  condition  and  perfonnance  assessment. 
As  the  leakage  current  increase,  the  size  of  the  dry  band  also  will  increase  leads 
to  the  dry  band  arcing,  deteriorate  the  insulator  perfonnance  and  contribute  to 
the  development  of  insulator  flashover.  Based  on  the  literature  study,  other  than 
the  existence  of  contaminations  on  the  surface  of  insulator  combined  with 
moisture,  the  variation  in  leakage  current  is  also  affected  by  the  environmental 
and  electrical  stresses.  Previous  researches  have  shown  the  effect  of 
environmental  and  electrical  stresses  on  surface  leakage  current  based  on 
experimental  and  simulation  results.  This  paper  outlines  an  analytical  approach 
based  on  dimensional  analysis  to  propose  a  new  mathematical  model  of 
leakage  current  under  environmental  and  electrical  stresses.  To  justify 
the  applicability  of  the  derived  dimensional  model,  the  new  model  has  been 
validated  using  previous  researcher’s  experimental  results.  The  validation 
indicated  that  the  proposed  model  had  shown  a  good  agreement  with 
the  previous  experimental  results.  The  proposed  dimensional  equation  for  this 
research  work  can  be  potentially  used  as  a  predictive  performance  model  to 
evaluate  and  monitor  the  leakage  current  and  insulator’s  perfonnance. 
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1.  INTRODUCTION 

Transmission  line  insulator,  also  known  as  overhead  insulator  is  greatly  affected  by  the  environmental 

stresses  produced  by  the  weather  and  surroundings.  Insulator  is  usually  exposed  to  uncertain  weather 

and  environmental  pollution  [1]  such  as  inert  material,  soluble  salt,  water,  electronic  conductive  dust, 
chemicals,  industrial  contaminants,  and  agricultural  pollutants  [2,  3].  Due  to  this  variety  of  pollutions, 
insulator’s  external  surface  becomes  contaminated  and  tends  to  accumulate  the  contamination,  allowing 
the  leakage  current  (LC)  to  flow.  The  contamination  deposited  on  the  insulator  surface  combined  with 
the  presence  of  moisture  dissolves  the  salt  deposited  and  encourages  the  formation  of  the  conductive  layer 
that  acts  as  a  highly  variable  and  nonlinear  resistor  to  be  formed  for  LC  to  flow  [3,  4].  With  high  LC  flowing 
between  the  high  voltage  and  grounded  side  of  insulator,  flashover  occurrence  that  will  lead  to  power  outage 
from  the  system  is  more  likely  to  happen  [5-7].  Flashover  phenomenon  reduces  the  power  system  reliability 
and  efficiency,  contributes  to  economic  losses,  and  permanent  insulator  impairment  [8,  9].  Due  to  these 
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reasons,  LC  monitoring  is  one  of  convincing  way  to  examine  the  insulator  performance  especially  under 
contaminated  conditions  [2,  10,  11].  LC  monitoring  describes  the  insulator’s  surface  condition  and  is  capable 
to  act  as  an  accurate  indication  tool  for  flashover  warning  and  hence  prolonging  the  insulator’s  lifetime  [12,  13]. 

In  recent  years,  researchers  have  investigated  the  effect  of  environmental  and  electrical  factors  on 
LC.  However,  most  of  the  research  results  come  from  experiments  and  simulations.  Testing  on  high  voltage 
prototype  sometimes  causes  costly  mistakes,  requires  a  relatively  long  considerable  time,  and  there  are  cases 
when  the  analytic  expression  of  the  variables  is  not  available  or  inaccurately  known  [14,  15].  Some  research 
efforts  were  also  directed  towards  establishing  a  mathematical  correlation  and  modeling,  which  is  adapted 
from  empirical  approach  [16-20].  However,  the  models  developed  only  effective  as  the  LC  prediction 
and  unable  to  represent  the  actual  physical  equation  that  explains  the  characteristics  of  the  LC.  This  paper 
presents  the  mathematical  relationship  to  describe  the  LC,  considering  the  environmental  and  electrical 
factors  using  the  analytical  approach  of  dimensional  analysis  (DA).  DA  is  an  analytical  technique  which 
applies  the  concept  of  similarity  to  explain  the  relationship  among  the  physical  quantities  [14,  15,  21]. 
By  using  DA  approach,  the  physical  phenomenon  is  represented  by  dimensionally  correct  equation  among 
certain  variables  [22,  23]. 


2.  DEVELOPMENT  OF  DIMENSIONAL  MODEL 

From  the  literature  studies,  it  can  be  concluded  that  the  LC  characteristics  have  a  strong  relation 
with  many  physical  variables  from  the  environmental  and  electrical  stresses.  The  LC  that  flows  on  insulator’s 
surface  depends  on  the  dominant  variables  such  as  Equivalent  Salt  Deposits  Density  (ESDD),  humidity, 
electrolyte  conductivity,  environmental  temperature,  applied  voltage,  environmental  pressure,  ultraviolet 
radiation,  and  insulator’s  creepage  distance.  A  mathematical  relationship  to  describe  the  LC  is  established 
between  these  variables  using  DA.  The  base  dimension  of  the  target  variables  is  shown  in  Table  1.  Among 
the  target  variables,  1L  is  considered  as  a  dependent  variable  and  ESD,  h,  aE,  Ta,  Vs,  pe,  Suv  and  Lc  as 
independent  variables.  The  unit  system  used  in  developing  the  LC  mathematical  model  is  SI  unit  standard 
system  and  the  base  dimensions  are  mass  (M),  length  (Z),  time  (7),  charge  (Q)  and  temperature  (0). 


Table  1.  Base  dimension  of  the  target  variables 


Variables 

Symbol 

Base  dimension 

Leakage  current  (LC) 

k 

QT-1 

Equivalent  salt  deposits  density  (ESDD) 

Esd 

MLT2 

Humidity 

h 

ML~ 3 

Electrolyte  conductivity 

°E 

M^L^TQ2 

Environmental  temperature 

Ta 

e 

Applied  voltage 

vs 

ML2T~2Q 

Environmental  pressure 

Pe 

ML -XT-2 

Ultraviolet  radiation 

Suv 

MT~ 3 

Insulator  creepage  distance 

Lr 

m~1l~1t2q 

From  the  target  variables,  the  relationship  between  the  independent  and  dependent  variables  are  established 
as  follows,  where  /  is  an  unknown  function. 

4  =  f  (ESD,h,aE,Ta,Vs,pe,  S-uv.Lc)  (1) 

From  the  list  of  the  target  variables,  the  number  of  variables,  Nv  =  9  and  the  number  of  dimensions, 
Nd  —  5.  DA  adopts  the  Pi  Buckingham’s  theorem  in  its  analysis.  Buckingham  method  starts  by  examining 
the  number  of  variables  involved  and  the  number  of  dimensions  [15].  The  final  result  is  a  relation  between 
Nv  and  Nd  which  provides  the  number  of  dimensionless  variables  (no  of  n’s).  Therefore,  by  applying  Pi 
Buckingham’s  Theorem,  the  no  of  independent  set  of  product  (n’s)  is  n’s=  Nv  —  Nd  =  9  —  5  =  4.  The 
dimensionless  product  can  be  written  in  a  form  as  in  (2). 

II  =  ESD£\h£2.aE£*.Ta£\Vs£\pe£*.Suv£\Lce*  (2) 

where  £± . £g  are  the  exponent  of  the  corresponding  target  variables.  The  dimensional  matrix  of 

the  variables  correspond  to  their  base  dimensions  with  MLTQ9  combination  is  arranged  as  follow; 
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The  rank  of  the  dimensional  matrix  is  established  as  5.  To  compute  the  independent  set  of  product, 
the  dimensional  matrix  is  partitioned  into  two  submatrices,  Matrix  A  and  Matrix  B.  Matrix  A  must  be  a  square 
matrix  and  not  singular.  If  Matrix  A  is  singular,  there  is  no  possibility  to  produce  the  independent  set  of 
products  (IT’s).  The  number  of  independent  set  of  product  (IT’s)  is  generated  by  using  the  formula  (3)-(4)  [15]. 


P  =  E.Z 


Z  = 


and 


E  = 


—A~1B 


where:  P 
B 
A 

q 

i 

6 

0 


:  Independent  set  of  product 
:  Matrix  B 
:  Matrix  A 

:  Exponent  of  dimension 
:  Identity  matrix 
:  Exponent  of  variable 
:  Null  matrix 


(3) 

(4) 


Matrix  P  established  in  (5)  described  the  target  variable  and  product  of  variables  [15]. 

1  0  0  0  " 

0  10  0 

0  0  10 

0  0  0  1 

0  0  0  0 


The  last  step  of  the  DA  modeling  is  to  produce  the  set  of  dimensionless  target  variables.  The 
complete  set  of  dimensionless  products  is  now  rewritten  as  in  (6)-(9). 


2 


„  _  >L-Pe5 

1*1—1  4 

(6) 

-Suv-L^ 

l 

( eSd)-suv2-lc 

U2=  1  „ 

(7) 

Vfpe  5 

h.Suv 2 

713  =  Pe3 

(8) 

1 

aEvs 

™  4  =  2  11 

(9) 

Pe5-^uv-^c 


By  referring  to  the  Buckingham’s  Theorem,  the  established  dimensionless  product  is  given  as  in  (10). 

nx  =  /(n2,n3,n3,n4)  (io) 

To  show  the  relationship  between  the  LC  and  other  independent  variables,  monomial  power 
form  [24]  is  applied  as  shown  in  (11).  Factors  a1,  a2  and  a3  is  assigned  by  considering  the  relationship 
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between  the  dominant  variables  and  leakage  current.  It  is  found  that  the  leakage  current  is  directly 
proportional  to  ESDD,  relative  humidity,  and  applied  voltage. 

7i i  =  /(7r2al.  7T3“2. 7r4“3)  where  aq  =  l,a2  =  1  and  a3  =  1  (11) 

Equation  (11)  is  rearranged  to  show  the  relationship  ofLC  to  the  other  dominant  variables. 

1 

h  =  P'5 f2Lcy  -f(n2al -n3a2 .7T4“3)  (12) 

Finally,  the  resulting  model  of  LC  characteristics  can  be  expressed  as  a  function  of  environmental 
physical  variables  in  (13)  where  Dc  is  the  dimensional  constant  and  can  be  defined  from  the  experiment. 
The  dimensional  constant  established  in  this  research  will  vary  depends  on  experimental  conditions. 

1 

h  =  Dc(Esd.  h.  aE.Suv4)  (13) 

From  (13),  it  can  be  proved  that  the  derived  dimensional  equation  obey  the  dimensional  homogeneity 
principle.  The  principle  of  dimensional  homogeneity  helps  in  describing  the  consistency  and  completeness  of 
an  equation  in  physical  algebra  [21].  Therefore,  (13)  is  projected  to  be  used  as  a  predictive  model  in 
estimating  the  LC  in  this  research. 


3.  VALIDATION  OF  THE  PROPOSED  DIMENSIONAL  MODEL 

For  validation  purposes,  the  proposed  leakage  current  dimensional  model  is  verified  with  other 
researcher’s  experimental  results.  The  value  of  the  dimensional  constant  or  test  parameter  is  recommended  to 
fit  the  proposed  model  with  other  researcher’s  experimental  results  by  using  Mathcad  Software.  This  section 
discusses  the  verification  of  the  model  at  different  ESDD,  electrolyte  conductivity,  relative  humidity,  applied 
voltage  and  insulator’s  creepage  distance. 

3.1.  Validation  of  the  model  at  different  electrolyte  conductivity 

Validation  of  the  proposed  model  at  different  electrolyte  conductivity  is  achieved  by  comparing  it  to 
the  experimental  results,  Suwarno  et  al.  [25]  on  polymeric  insulator  subjected  to  inclined  plane  test  under 
various  levels  of  kaolin-salt  pollutions.  All  the  dominants  variables  (ESD,h,Vs,  pe,Suv,  Lc)  in  (13)  are 
assumed  to  be  constant  except  the  levels  of  electrolyte  conductivity.  By  assuming  all  the  mentioned  variables 
are  constant,  the  proposed  model  can  be  simplified  as  follows; 

IL(a)  =  D1.a  +  I0  (14) 

where  D±  is  the  recommended  dimensional  constant  to  fit  in  the  model  with  experimental  results.  Equation  (14) 
is  altered  by  adding  I0  to  adapt  the  experimental  result.  Based  on  the  simplified  model,  it  is  found  that 
the  LC  flows  on  insulator  surface  is  significantly  increase  with  electrolyte  conductivity.  Figure  1  shows 
the  comparison  of  the  LC  from  the  model  developed  and  the  doted  data  from  the  experimental  result  [25]. 
It  can  be  observed  that  the  proposed  model  fitting  is  very  close  to  the  experimental  results.  Table  2  shows 
the  proposed  LC  equation  in  relation  to  the  electrolyte  conductivity  with  its  data  error  and  standard  deviation. 


Electrolyte  Conductivity  (mS/cm) 


Figure  1.  Comparison  ofLC  of  experimental 
and  proposed  model  at  different  electrolyte 
conductivity  [25] 


Table  2.  Empirical  equation  ofLC  in  relation  with 
electrolyte  conductivity 


Empirical  equation 

=  0.72(7,;  +  1.2 

Error 

0.0592 

Standard  deviation  (Model) 

1.0824 

Standard  deviation  (Exp) 

1.0601 
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From  Table  2,  it  shows  that  the  model  has  a  quite  low  deviation  of  5.92%  with  reference  to 
experimental  surface  LC.  By  varying  Du  the  curve  can  be  adjusted  to  suit  the  experimental  data.  This 
provides  an  approximation  to  form  the  equation  between  LC  and  ESDD  as  in  Table  2. 

3.2.  Validation  of  the  model  at  different  ESDD  and  applied  voltage 

The  proposed  model  is  verified  with  the  experimental  data,  Banik  et  al.  [26],  at  different  ESDD 
and  applied  voltage.  Experiment  was  performed  on  porcelain  insulator  to  measure  surface  LC  with  different 
pollutions  conductivity  and  increasing  of  applied  voltage  in  steps  of  5  kV  to  maximum  40  kV  under  control 
environmental  [26].  The  independent  physical  variables  except  for  ESDD  and  applied  voltage  are  assumed  to 
be  constant.  Hence,  (13)  can  be  formulated  in  a  function  of  ESDD  and  applied  voltage. 

Il(Esd,Vs)  =  D2.Esd.VsI  +  I0  (15) 

By  referring  Figure  2,  it  can  be  observed  that  the  model  is  in  complete  agreement  with  the  previous 
experimental  results  [26]. 


| 

I  ila 

B  *  *  * 

V  IL2(ESDD) 


(a) 


(b) 


Figure  2.  Comparison  between  model  and  experimental  result  of  LC  at  different  ESDD 
and  applied  voltage  [26],  (a)  10  kV,  (b)  20  kV 


The  model  able  to  compute  the  LC  for  different  ESDD  and  applied  voltage  with  deviation  of  6.04% 
and  11.86%  for  10  kV  and  20  kV  of  applied  voltage  respectively  (refer  Table  3).  The  LC  is  in  linear 
relationship  with  ESDD  and  applied  voltage.  Increasing  the  applied  voltage  will  increase  the  electric  field 
distribution  and  affect  the  LC  variation  [27,  28].  The  dimensional  constant  of  D2  also  shows  an  increment  in 
its  value  with  increasing  of  the  applied  voltage. 


Table  3.  Empirical  equation  of  LC  in  relation  with  ESDD  and  applied  voltage 


Empirical  equation 

h  =  0.57 Esd.  Vs 1  +  0.18 
(Applied  voltage  10  kV) 

/l  =  3.15(£ot).KsI  +  0.28 
(Applied  voltage  20  kV) 

Error 

0.0604 

0.1186 

Standard  deviation  (Model) 

0.0956 

0.6066 

Standard  deviation  (Exp) 

0.1118 

0.7127 

3.3.  Validation  of  the  model  at  different  relative  humidity  and  ESDD 

Equation  (13)  is  written  in  a  function  of  relative  humidity  and  ESDD  as  shown  below.  From  (16),  it 
is  found  that  LC  is  in  linear  proportion  with  relative  humidity  and  ESDD.  LC  increases  with  the  increasing  of 
contamination  severity  and  relative  humidity.  Increasing  of  contamination  severity  combine  with  high 
relative  humidity  add  up  the  moisture  content,  reduce  the  insulator  surface  resistance  and  lead  to  a  high  LC 
flows  on  insulator  surface  [29]. 

IL(h,  Esd)  =  D3 .  Esd ,h  +  I0  (16) 

The  proposed  model  is  verified  with  experimental  results,  Tousi  et  al.  [2],  Figure  3  shows 
the  verification  of  the  proposed  model  using  the  data  from  the  experiment  [2],  It  can  be  observed  that 
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the  deviation  of  the  model  is  quite  low  (see  Table  4)  which  are  4.51%  and  4.01%  compared  to  the  experimental 
result.  This  identify  that  the  model  is  capable  to  predict  the  LC  at  different  relative  humidity  and  ESDD. 


Figure  3.  Verification  of  proposed  model  at  different  relative  humidity  and  ESDD  using  the  results  from 
experimental  [2],  (a)  ESDD  0.26  mg/cm2,  (b)  ESDD  0.358  mg/cm2 


Table  4.  Empirical  equation  of  LC  in  relation  with  relative  humidity  and  ESDD 


Empirical  equation 

IL  =  1.71Esd.  h  +  0.12 
ESDD  0.26  mg/cm2 

4  =  1.364sD.  h  +  0.16 
ESDD  0.358  mg/cm2 

Error 

0.0451 

0.0401 

Standard  deviation  (Model) 

0.0889 

0.0974 

Standard  deviation  (Exp) 

0.1102 

0.1124 

3.4.  Validation  of  the  model  at  different  insulator’s  creepage  distance 

To  validate  the  LC  model  at  different  insulator’s  creepage  distance,  (13)  is  simplified  in  a  function 
of  specific  creepage  distance  as  shown  in  (17)  by  assuming  all  independent  variables  except  creepage 
distance  as  constant. 

6 

/l(Lc)  =  D4.l7  +  /0  (17) 

The  model  is  compared  to  previous  research  data,  Parihar  et  al.  [30],  which  presents  the  leakage 
current  on  composite  long  rod  insulator  with  different  specific  creepage  distance  and  contamination  severity. 
Figure  4  shows  the  LC  in  dependency  on  the  specific  creepage  distance  of  the  proposed  model  and 
the  previous  research  data  [30].  From  the  graph,  it  can  be  observed  that  the  LC  is  inversely  proportion  to 
the  specific  creepage  distance.  With  a  higher  value  of  specific  creepage  distance,  conducting  layer  could  be 
prevented  thus  decreases  the  LC  flows  on  insulator  surface.  From  Table  5,  it  may  simply  have  verified  that 
the  model  has  a  very  close  agreement  with  the  previous  experimental  results  [30]  with  a  deviation  of  0.23%. 


Specific  Creepage  Distance  (min/  kV) 

Figure  4.  Comparison  of  model  and  previous  research 
result  at  different  insulator’s  creepage  distance  [30] 


Table  5.  Empirical  equation  of  LC  in  relation  with 


insulator’s  creepage  distance 


Empirical  equation 

6 

4  =  147.3L/  +  4.2 

Error 

0.0023 

Standard  deviation  (Model) 

0.3783 

Standard  deviation  (Exp) 

0.3702 
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4.  CONCLUSION 

The  present  article  describes  an  analytical  approach  on  modeling  the  LC  of  insulator  under 
environmental  and  electrical  stresses  based  on  DA  technique.  The  proposed  dimensional  model  derived  in 
this  article  complies  with  homogeneity  condition  and  expressible  by  a  dimensionally  homogeneous  equation 
in  terms  of  specified  parameters.  The  proposed  method  also  able  to  present  the  dominant  parameters 
associated  to  the  LC  behavior  and  recommends  the  appropriate  correlation  of  the  test  parameter  value  to  be 
used  when  executing  the  test.  Throughout  this  paper,  it  has  been  observed  that  the  dimensional  model  has 
a  good  agreement  with  the  experimental  results.  This  clearly  shows  that  the  model  able  to  estimate 
the  surface  leakage  current  on  insulators  under  the  influence  of  environmental  and  electrical  stresses 
and  contribute  to  an  indication  mechanism  of  flashover  warning.  Furthermore,  the  proposed  work  also  might 
help  in  reducing  the  cost  of  testing  and  labor  work  for  the  study  of  the  insulator  subjected  to  environmental 
and  electrical  factors. 
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